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Abstract 

We give a path integral expression for the quantum amplitude to produce a black 
hole from particle collisions. When expanded about an appropriate classical solution 
it yields the leading order contribution to the production amplitude in a curvature 
expansion. Classical solutions describing black hole production resulting from two par- 
ticle scattering at non-zero impact parameter, combined with our formalism, indicate 
a geometric cross section for the quantum process. In TeV gravity scenarios these so- 
lutions may exhibit large curvatures, but (modulo a mild assumption about quantum 
gravity) corrections to the semi-classical cross section are small. 
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1 Introduction 



Recently proposed models with extra dimensions solve the hierarchy problem by bringing 
the fundamental scale of gravity (henceforth referred to as the Planck scale) down to the 
electroweak scale |lj . In such scenarios, quantum gravitational effects arise at energies as low 
as a TeV. Perhaps the most dramatic example of such phenomena is black hole production 
in particle collisions with center of mass energy greater than the Planck scale 0-0- Such 
events would lead to dramatic signatures at colliders and in cosmic ray collisions, and perhaps 
imply an end to our ability to probe shorter and shorter length scales @]-[[{J. 

The cross section for black hole production in high energy collisions is difficult to compute. 
In [[J-0 it was asserted that the cross section is geometrical, determined by the impact 
parameter at which the particle pair at closest approach is within the Schwarschild radius 
associated with the center of mass energy ^/s. If this is the case, black holes would be 
copiously produced at LHC and in cosmic ray collisions. However, Voloshin has criticized 
these claims, arguing for an exponentially small cross section ||. 

Eardley and Giddings || have analyzed classical solutions in general relativity which 
describe two particle high energy collisions at non-zero impact parameter (see also |10|-|14||). 



They demonstrate the existence of a closed trapped surface for any collision with sufficiently 
small impact parameter (at fixed center of mass energy). Their lower bound on the critical 
impact parameter leads to a geometrical classical cross section in rough agreement with the 
earlier naive estimates. 

The solutions in || yield immediate answers to Voloshin's two main objections: (see also 
[|15] for a model calculation that shows no exponential suppression) 



1) Euclidean suppression: because there are classical trajectories with two particle initial 
conditions which evolve into black holes, the process is clearly not classically forbidden, and 
hence there is no tunnelling factor. 

2) CPT (time reversal): Voloshin argues that since black holes produce a thermal spec- 
trum of particles during evaporation, rather than a few highly energetic particles, the time- 
reversal of the production process (and hence the production process itself, by CPT) must 
have very low probability. However, the time-reversed classical solutions exhibit a very ener- 
getic wave of gravitational radiation colliding with the time-reversed black hole to produce 
the two particle state (in the formation process this is the energy which escapes the hole). 
The process is not thermal and involves very special initial (final) conditions. 

In this letter we use a path integral formalism previously developed by Gould, Hsu 
and Poppitz JR], [17j (GHP) to study non-perturbative scattering. In the GHP formalism 



2 



quantum S-matrix elements are computed in a systematic expansion about classical solutions 
satisfying appropriate boundary conditions. We show that the formalism is readily adapted 
to the problem of black hole production from high energy collisions. 

Although some black hole production amplitudes (for example, involving many soft initial 
particles) can be computed unambiguously from classical solutions, in the two particle case 
quantum corrections might be significant due to large curvatures. However, we argue that 
even in this case the semi-classical approximation is probably a good one. We explain the 
importance of large curvatures to the quantum corrections from an effective field theory 
point of view. 



2 Path Integral Formalism 

In this section we review a method formulated to describe scattering processes involving 
non-perturbative field configurations using the path integral [16, Previous interest in 
this problem centered on baryon number violation in the electroweak theory. In this context, 
it was realized that classical configurations satisfying vacuum boundary conditions (such as 



instantons) were unsuitable for the computation of high energy scattering A good 

semiclassical approximation in this context requires taking into account the initial state 
overlap as well as the classical action. 

In [0], an exact expression is derived for the S-matrix describing a transition from an 
initial two particle state to a coherent final state. (See below for a precise definition of co- 
herent state.) This expression is approximated in the usual stationary phase approximation, 
which leads to a boundary value problem involving the usual classical equations of motion, 
but with boundary conditions determined by the initial and final states. Unfortunately, nu- 
merical searches for "interesting" solutions of this boundary value problem in lattice gauge 
theories have only uncovered configurations describing transitions between states consisting 



of many soft quanta Hl9] . The Eardley-Giddings result is the first we are aware of in which 
a high momentum initial state evolves into a large, low momentum final state. However, we 
should note that the two particles in the initial state are "dressed" by strong gravitational 
fields, so the effective number of gravitons in the initial state is actually large. 

We first give a brief review of the general method, before addressing issues particular 
to general relativity and black holes. In the path integral representation of the process 
N) — *■ trajectories are weighted by the appropriate action exp^S], as well as by the 
overlap of the asymptotic part of the path with the initial and final states. Therefore we 
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expect that the amplitude must be expressible in the form |17|: 



(/i5ii> ~ j dfa d4 f d<p nm)} co 

where we have explicitly indicated the fluctuations of the fields in the asymptotic past and 
future (Tij) in the measure. The wave-functionals VPjj measure the overlap of the initial 
and final states with eigenstates at asymptotic times. In [16 a derivation of (U) was given, 



along with the explicit form of the wavefunctionals in the case that the initial and final states 
are either wave-packets (including plane waves) or coherent states. 

In the limit Tjj — > ±oo, the amplitude in (fj) is just the S-matrix, S/j. The GHP 
procedure amounts to a semi-classical evaluation of this object taking into account the initial 
and final state overlaps in the extremization. The result is a boundary value problem with 
boundary conditions determined by the initial and final state, but governed by the usual 
equations of motion. Here our goal is to calculate the S-matrix element between an initial 
two particle state and a final state which includes a black hole. 

First, we express the kernel of the S-matrix in a basis of coherent states (see, e.g., the 
text by Fadeev and Slavnov P0[). The initial and final states are defined by sets of complex 
variables a = {a^}, b* = {&£}, respectively. A coherent state |ctk) is an eigenstate of the 
annihilation operator a^. d^a^) = ^kl^k)- Recall that coherent states saturate minimum- 
uncertainty bounds, and hence are good (although not unique) candidates for semi-classical 
states. 

The transition amplitude from an initial coherent state |a) at time Tj to a final coherent 
state |b*) at time T/, can be expressed 

(b*|C/|a) = Jd&dtPj (b*\<p f ) ((PflUl^) i&la) , (2) 

where U is the evolution operator between time Tj and Tf. A "position" eigenstate of the 
field operator is denoted | 0) and <pij = <f>(Tij). Then, from @, we obtain the S-matrix 
kernel in a compact form in terms of path integrals 

(b*|£|a) = 5[b*,a] = lim f e% e% e B ? e B > [*' D<\> e iS ^ , (3) 

Tt,Tf -*^foc J Jfc 

where S [0] is the action functional. The path integral appearing here is over fields obeying 
the boundary conditions 4>{Tij) = (pij. The functional Bf is 

B f [b*,<p f ] = -\ Jd 3 k b* k b*_ k e 2 ^ r / - i Jd 3 k o; k 4f(fi M-k) (4) 
+ J d 3 k e iuJ * T f 6* <t> f (-%) , 
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in terms of which the wave functional of the final coherent state is ( b* | 0/ ) = exp (Bf [b*, </>/] 
Similarly, the functional Bi can be expressed in terms of the initial coherent state wavefunc- 
tional (0j|a) = exp ( Bi [a, ) , with o k replaced by and Tf by — Ti in (Q). The 
3-dimensional Fourier transform is defined 

4>i, f (k) = j^m^^j^)- (5) 

The kernel (|3]) is a generating functional for S-matrix elements between any initial and 
final iV particle states, by functional differentiation with respect to arbitrary and o k . We 
now use this fact to construct a kernel for scattering from initial two particle states. We 
define an initial two particle (wave packet) state at t = Tj 



p, -p 



) = J d 3 k a R (k) d[ J d 3 k' a L (k') a{, | 0) , (6) 



where a k is a creation operator, and a R) L{k) are arbitrary smearing functions of k, localized 
around some reference momenta p and — p respectively. The wave packets are normalized so 
that / d 3 k | a R>L (k) | 2 = 1 . 

This state can be generated by functional differentiation of the coherent state |a) with 
respect to ak 



\P, ~P) = [ d 3 k d 3 k! a R (k) a L (k') — - — | a) 
J oat oa-^i 



(7) 

a = 



So, differentiating under the functional integral, the S-matrix element between the two 
particle state (|7]) and any final state |{0kJ) involves the following functional at t = Ti 

5 5 , r 

exp ( Bi [a, (pi 



a=0 



5ak 5a^ 

2 V^J^ <pi(k) <Pi(k') e-^ +w ^ exp (5, [0 



) VI J J 1 



after dropping a term which vanishes in the limit Tj — > — oo. The last factor here is simply 
the normalization of the initial position eigenstate 

exp ( -- J d 3 k u k <pi{k) <pi{-k) \ . (9) 

We combine this with the smearing functions and finally obtain an S-matrix kernel for 
the scattering of two wave packets into arbitrary final states, 

S [b*,2] = lim / dfadfa a R -fa at L - <p t e B A b M + B ^ f* D<\> e lSW , (10) 
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where we have denoted the initial state ([F]) by "2" . Here we have used the following compact 
notation for the initial state factors 



a 



d 3 k a(k) fa(k) . (11) 



We exponentiate the initial state factors into an "effective action" , so that 

S[b*,2] = lim I ' d(f> f d(j>iD(j> e r , (12) 

where the effective action V is 

T [0] = lna R -faa L -fa + Bi[Q,fa] + iS[(j)\ + B f [b*,4> f ] , (13) 

after dropping a term which vanishes as Tj — > — oo. 

We can now derive the boundary value problem by varying the effective action. Varying 
the entire exponent T with respect to 0(x) for Tj < t < Tf gives the source-free equations of 
motion 

M . (14) 



5(f)(x) 

Varying the entire exponent with respect to <pi{k), gives 



ifa(k) + u k Hk) = V2^ ( + ^#1 e-*** 



ocr ■ (pi a L - 



(15) 



The first term on the left hand side comes from a surface term in the action S. The other 
terms come from variation of the wave functional at t — Tj. This boundary condition involves 
both the positive and negative frequency parts of the field. 



The boundary condition fll5|) at the initial time slice is rather complicated. However, it 
can be simplified since a real field may be written in the asymptotic region t = Tj —>■ — oo 
clS db plane wave superposition 

1 



fa(k) = -= (uue-^ Ti + u*_ k e^ Ti ) . (16) 



Equation (O) then reduces to the requirement 



a R (k) + a L (k) 

M k = T\T72 ' ( 1T ) 

1 + / d 3 k aR^k) ai(k) 



using the normalization condition on oc^r. This solution is consistent with physical intuition, 
the classical field reducing to the initial particles at early times. The overlap of the left- and 
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right- moving wave packets in the denominator is very small for narrow high energy wave 
packets. 

A similar analysis relates the late time boundary condition on <p to the coherent state 



b* [16]. The classical field satisfying these boundary conditions extremizes the S matrix 
for production of the coherent state in a two particle collision: 5[b*,2]. Note that for 
arbitrary choice of b* there is no guarantee of a classical solution satisfying the necessary 
boundary conditions: in some cases a complex trajectory extremizes the S-matrix, leading 
to exponential suppression of the process |pL6| . However, conversely every classical solution 



obeying initial conditions (|T6|)-(|T7j) corresponds to an unsuppressed quantum amplitude. 

Now consider the extension of this formalism to general relativity, and to black hole 
production. Clearly one can replace <ft with the metric field g^ u plus appropriate matter 
fields. While general covariance does not permit a unique time slicing, the gravitational 
action S g = J d^Xy/^g R is still well-defined and, in fact, the action can be expressed in 
terms of a surface integral over Bondi masses ET1 . 



The notion of an S-matrix is appropriate if we consider asymptotically flat spacetimes in 
the far past and future, plus additional excitations. Black holes themselves are considered 
excitations, and we must extend our Hilbert space to include quantum states representing 
black holes. A pragmatic way to approach this is to define (semi-classical) black hole states 
as those with a strong overlap with the trajectories corresponding to classical black holes. 
(These must of course have mass much larger than the Planck mass.) In a classical black hole 
solution, excess energy is radiated away by late times and the exterior metric can be classified 
by a limited number of quantum numbers such as mass, charge and angular momentum. A 
minimal formulation involves a Hilbert space of black holes classified by their exterior metric 
at future null infinity. 

Although it is beyond the scope of this letter, it is worth noting that a more detailed 
analysis of black hole states, which takes account of the internal structure of the black hole 
(i.e. the fields inside the horizon, not just at future null infinity), indicates more states 
than are counted by externally visible quantum numbers. That is, if one considers the 
internal structure of the black hole in defining the Hilbert space, there are many additional 
semi-classical states that one might associate with a given exterior metric, but which differ 
radically within. Gedanken experiments involving semi-classical black holes in this formalism 
might teach us something about black hole information. For example, relative phases and 
interference patterns due to internal structure might be observable in black hole scattering. 

Any classical solution connecting particle-like initial conditions to an asymptotic black 
hole configuration provides an extremal configuration about which to expand the S-matrix; 
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the leading contribution is a pure phase with no exponential suppression |T(|. In the next 
section we consider quantum corrections to this leading semi-classical approximation. 



3 Quantum corrections 



In the original application of |L6| to quantum fields in flat spacetime, it was shown that the 
quantum corrections to the semi-classical approximation to the S-matrix are suppressed by 
powers of the coupling constant. In particular, the corrections can be expressed in terms of 
propagators and interaction terms which result from expanding about the classical solution. 
All interaction terms carry explicit powers of the coupling, resulting in a well-defined loop 
expansion. 

In general relativity there is of course no small dimensionless parameter. Expanding 
about a background configuration yields interactions which are suppressed by the background 
curvature in Planck units. Classical solutions describing the ordinary gravitational collapse 
of many "soft" particles (e.g., collapse of a large star or dust ball) can produce black holes 
without regions of large curvature. Our formalism applies directly to such solutions, resulting 
in a semi-classical amplitude without large quantum corrections. 

In the two-particle solutions of |§, regions of large curvature can arise quite early in 
the evolution (e.g., when shock fronts collide), even if the black hole produced is large 
(y/s » Mpi anc k). If one takes the size of the colliding particles to be of order the Planck 
length L, one finds curvatures at the shock front of order s. In this case, quantum corrections 
might be large. In fact, we run into a fundamental problem concerning quantum gravity. 
Because gravity is non-renormalizable, we have to consider all possible generally covari- 
ant higher dimension operators in our lagrangian, such as higher powers of the curvature. 
Certainly, such terms will arise from the ultraviolet part of any loop calculation and will 
presumably be only partially cancelled by counterterms. In an effective lagrangian descrip- 
tion, we expect these operators to be present, but suppressed by powers of the Planck scale. 
In large curvature backgrounds these terms may not be negligible, so the size of quantum 
corrections will in principle depend on unknown details of quantum gravity. 

We can state this conclusion in a slightly different way. Consider the classical evolution 
of some initial data in an effective low-energy description of gravity. We can only trust the 
Einstein equations (which result from the lowest dimension term in the effective lagrangian) if 
large curvatures are never encountered during the evolution. Once a region of large curvature 
is encountered, subsequent evolution might depend in detail on the nature of the higher 
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dimension operators, and hence on the nature of quantum gravityQ In it is suggested that 
since a closed trapped surface is identifiable in a low-curvature region, the classical solution 
is a good guide to the true quantum behavior. Strictly speaking, this is not sufficient - the 
solution has already evolved through a potentially high curvature shockwave region. 

It is reasonable to expect that short distance features of the metric which lead to high 
curvatures do not affect the classical Einstein evolution of the solutions in || on large 
length scales, such as of order b, the impact parameter. The large distance behavior of 



the Aichelburg-Sexl metric |LC| used in U is independent of the short distance features of 
the particle as long as its size r is much smaller than b. In fact, Kohlprath and Veneziano 
have recently extended the Eardley-Giddings construction to the case of colliding particles 
of finite size 1)231. 



In TeV gravity scenarios r is likely to be the Planck length. However, we can instead 
consider collisions of particles of size much larger than the Planck length but much less than b. 
(In other words, "colliding Jupiters" at ultra-relativistic velocities, with b much larger than 
the radius of Jupiter!) By adjusting the impact parameter and particle size relative to the 
Planck length (while keeping a large hierarchy between the three) we can keep the shockwave 
curvatures parametrically smaller than M| lanck while preserving the long distance behavior 
that leads to horizon formation. In Planck units, the maximum curvature in the shockwave 
is R ~ (6/r) 2 (L/r) 2 , where h ~ L 2 yfs and L is the Planck length. By independently varying 
-y/i and r, we can make the first ratio large and the second small, while keeping their product 
small. The semi-classical approximation applies quite well in this limit and higher dimension 
operators can be neglected. 

Barring unexpected quantum gravitational effects which are sensitive to the size of the 
objects colliding (i.e. that make long distance behavior dependent on the short distance 
metric), the quantum cross section for black hole production will be well approximated by 
the semi-classical one even in TeV gravity scenarios. 
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